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SUMMARY 


A procedure is presented in the form of charts which permits the 
rapid estimation of the natural bending frequencies of helicopter rotor 
blades, both rotating and nonrotating. Since the approach is based on 
Southwell’s equation,’ an evaluation of the method with regard to such 
things as higher modes, blade offset, and variable mass and stiffness 
distributions is also given. The evaluation shows that, when nonrota- 
ting beam bending modes are used, Southwell's equation yields reason- 
ably accurate bending frequencies for rotating helicopter blades. 

Example conrparisons of frequencies estimated using the charts with values 
given by the manufacturer for several actual blades show that the sim- 
plified procedure yields good practical results . 


INTRODUCTION 


The purpose of this paper is to present results in chart form which 
pennit the rapid estimation of bending frequencies of rotor blades . The 
proposed method of frequency determination maices use of the familiar 
Southwell form; thus, it is also the purpose of this paper to show that 
this approach works quite well when such things as higher modes, blade 
offset, and variable mass and stiffness distributions are considered. 

The paper is divided into three parts as follows ; In the first part a 
review and an evaluation of the Southwell approach is given, in the sec- 
ond part frequency charts are presented, and in the third part the 
results of applying these charts to some actual helicopter blades are 
given to Indicate the order of accuracy obtainable in practical cases 
with this procedure. 


REVIEV7 AND EVALUATION OF SOUTHWELL APPROACH 


The basic form of the Southwell eqmtion which defines the bending 
frequencies of a rotating beam can be obtained by energy considerations 
and is given in reference 1 as 
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where 

0 ^ bending frequency of rotating beam 

ft rotational speed of beam 

m mass distribution for beam 

El stiffness distribution for beam 


T]_ft^ tension force in beam 

L length of beam 

X spanwise coordinate along beam 

y>y'#y" beam mode shape and derivatives with respect to x 

Eqtiation (l) yields exact values for the bending frequencies of a 
rotating beam if the mode shapes of the rotating beam are known. Since 
these exact shapes are usually not known, however, it is necessary to 
assimae mode shapes from which approximate frequencies can then be esti- 
mated. If the nonrotating mode shape is substituted into the first 
term of equation (l) and the coefficient of ft^ is replaced by K, the 
Southwell constant, the Southwell eqixation takes the following form: 


Kft^ 


( 2 ) 


where 

cofjp nonrotating bending frequency of beam 

K Southwell constant 

If the hinge or point of fixity of the rotating beam is offset from the 
axis of rotation, it may easily be shown that the Southwell constant can 
be written in the form: 



NACA EM L54G02 


5 


K = Kq + Kj_e 


(5) 


where 

Kq zero -offset Southwell constant 

offset correction coefficient for Southwell constant 

e distance hinge is offset from axis of rotation, percent 

beam length 

If equation (5) is substituted into equation (2), 

0^2 = + (Kq + Kie) (^4-) 

If certain constants for th.e 1368111 are introduced into equation (4), it 
can be written as 


( O ^ = + ^Kq + Ktlc) (5) 

moL^ 

where Iq and idq are measured at the root of the beam and 

L length of beam outboard of hinge or point of fixity 

a^i nonrotating frequency coefficient for beam vibrating in 

nth mode 

In order to provide a basis for estimating the accuracy, usefulness, 
and possible limitations of the Southwell approach, a systematic series 
of blades was selected and frequencies were calculated by the Southwell 
approach and a more exact process. Figure 1 shows the cases studied using 
both methods . These cases incliBie the uniform cantilever with O, 5^ 

10 percent offset and the viniform hinged beam with the same offsets . The 
"linear" tj-pe beams are beams in \dilch the mass and stiffness both vary 
linearly from the root value to zero at the tip. Zero- and 10-percent 
offset were treated for this type of beam with both the cantilever and 
hinged root . 

The Southwell frequencies were obtained with the iise of the inn <^ e 
shapes for the nonrotating beam whereas the true or reference frequen- 
cies were obtained by a Rayleigh-Ritz energy procedure involving 
expansion of the rotating beam modes in terms of the nonrotating modes 
of a uniform beam. For the cantilever beams, five nonrotating uniform 
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nodes were used in the expansion. For the hinged beams, a pendulm mode 
was used in addition to xhe five nonrotating nodes . Some of the results 
of this investigation are shown in figures 2 to 5* 

The variation of bending frequency with rotational speed for a 
uniform hinged beam is shown in figure 2. The range of 

in this figure is rouglily 25 percent above that encountered in current 
helicopters . From the resiilts shown it is evident that the Southwell 
results are quite accurate, the max:imimi error being about 5 percent in 
frequency sqxiared or only about 1^ percent in frequency. This maximum 

error is about the sane for all three modes. 

In order to avoid confusion, res\ilts for the 5-pei’cent-offset case 
are not sha^m in this figure. They fall roughly midway between the 
0- and 10-percent -offset curves and show the same type of agreement 
between exact and Southwell results . 

Frequency results for the "linear" type hinged beam are shown in 
figure 5 * From this figure it is apparent that the Southwell results 
are extremely acc\rrate, even for the highest rotational speeds shown. 

A comparison of frequency results for the uniform and the "linear" 
type hinged beam is given in figure 4 . The most important thing to be 
noted from this comparison is the difference in slope between the dashed 
and solid line for each mode. The slope of each of these lines is 
directly proportional to the Southwell constant. The large difference 
in slope, particularly evident for the first mode, indicates that a 
single value of the Southwell constant for each mode could not ade- 
quately predict the varlaxions of frequency with rotational speed which 
are shown in the figure. 
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1st 


Frequency results for uniform cantilever beams are presented in 


figure 5 • The range of 


-|2 


^ist 


corresponds roughly to that covered 


for the hinged beams ; the difference in scale is due to the fact that 
the first bending frequency of a uniform hinged beam is about four times 
greater than the first bending frequency of a uniform cantilever beam. 
For each mode the lower dashed and solid curve are for zero offset and 
the upper pair are for 10-percent offset . 


From the figure it is apparent thaL the Southwell results are very 
accurate for the second and third modes. For the first mode, however, 
the error in frequency is somewhat larger, about 5 percent. Similar 
resiilts obtained for cantilever beams with mass and stiffness distribution 
varying linearly from the root value to zero at the tip show about the 
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sane type of agreement Taetween the Southwell and the more exact results 
as are given in figure 5* Since the error in Southwell results for the 
first mode is roiighly the same for hoth these cases it is reasonable to 
assume that applying a correction factor determined by the error shown 
in figure 5 will result in more accurate prediction of the first mode 
bending frequency of other cantilever beams as well. 

From the foregoirsg evalviation of the Southwell approach it was con- 
cluded that Southwell constants based on the nonrotating beam mode shapes 
lead to reasonably accurate bending frequencies of rotating helicopter 
blades. The evaluation also showed that the Southwell constants vary 
appreciably with beam mass and stiffness distribution. 


CEAETS FOE FREQUENCY DETERMINATION 


In order to provide a meems for rapidly estimating rotor-blade 
bending frequencies, the nonrotating frequency coefficients, zero-offset 
Southwell constants, and offset correction coefficients for Southwell 
constants have been computed for a series of beams with linear mass and 
stiffness distributions. The range of mass and stiffness distributions 
was selected to encompass variations found in currently manufactured 
blades with some latitude for new design. All the constants are based 
on the mode shapes of the nonrotating beam which were obtained by stand- 
ard numerical iteration procedures . The cantilever modes were obtained 
using 10 stations while the hinged modes were obtained using 15 stations. 


The variation of the nonrotating frequency coefficient a^, with 

beam mass and stiffness distribution is shown in figure 6. The abscissa 
is the ratio of tip mass to root mass, 1.0 represents a constant mass 
beam, and 0 the case where the mass varies linearly to zero at the tip. 

The solid curves are for beams with constant stiffness along the 
length, the dashed curves for beams where the stiffness drops to half 
the root value at the tip, and the long- and short -dashed curves for 
beams where the stiffness is zero at the tip. 

Figirre 7 permits selecting an accompanying value for the Southwell 
constant for beams with zero hinge offset. The variation with stiffness 
distribution appears quite small, particularly for the first mode. 
Actually, however, there is a maximum difference of about 5 percent. 

The variation with mass distribution is obviously somewhat more pro- 
noimced for all the modes shown. The zero -offset Southwell constant is 
not shown for the zero or pendulum mode since it is always unity inde- 
pendent of the mass and stiffness distribution of the beam. 
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In order to account for cases where the hinge is offset from the 
axis of rotation, offset correction coefficients Kq_ have also been 
computed for this family of beams. These coefficients, when miiltiplied 
by the offset given as a fraction of the free beam length, yield the 
correction to be added to the zero-offset Southwell constant. 

The variation of this offset coefficient is shown in figure 8. It 
is evident that this variation is quite similar to that shown in figure 7 
except for the zero mode. An expression for the zero-mode offset coef- 
ficients presented here is given in reference 2. 

Results similar to those presented in figures 6 to 8 have also been 
obtained for cantilever beams and are shown in figures 9 to H- The 
effect of root fixity on the Southwell constant for the various modes can 
be deduced by comparing the curves in figures 8 and 9 with those in 
figures 10 and 11. It should be mentioned that a more accurate estimation 
for the first bending frequency of a rotating cantilever beam should be 
obtainable by reducing the frequency obtained usings the charts and Southwell's 
equation by a small percentage which may be quickly estimated from figure 5* 


EXAMPLE RESULTS 


In order to Illustrate the type of accuracy which can be expected 
in using the frequency charts of figures 6 to 11, bending frequencies 
have been estimated for the first three modes of four existing helicop- 
ter .blades, all of which are hinged. The following procedure was used 
in this estimation: 

(1) Straight lines were faired through the m and El distribu- 
tions for the blade, large vsQ.ues neax the root being ignored. 

(2) From these fairings, the effective root values of m and El and 
the necessary tip-to-root ratios were obtained. 

(5) By using these ratios, values of a^^, Kq, and K]_ were obtained 
from the charts . 

(4) Substitution of these constants into the Southwell equation 
yielded the bending frequencies at zero and at rated rotor speed. 

The results are shown in table I. The m and El distributions 
for the blades are shown on the left-hand side of the table. The actual 
distribution is given by the solid lines. The dashed line is the linear 
approximation selected to represent this variation. It sho\ild be empha- 
sized that these linear approximations used in estimating the freqtien- 
cies were the initial ones selected and were not juggled to obtain the 
best agreement. The frequencies shown as exact in this figure are values 
furnished by the manufacturer. 
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If the exact and estimated results for the "blades are compared, it 
is evident that the resiilts are quite accurate when the crudeness of 
the linear approximations used is considered. It is interesting to note 
that in all cases the estimated frequency of the rotating beam is more 
accurate than that for the nonrotating beam; this indicates that the 
linear approximations yield more accurate values for the Southwell con- 
stants than for the nonrotating frequency coefficients . 

Although no comparisons have been made for fixed end blades, it is 
believed that even more accurate results should be obtainable for this 
end condition since large values of root stiffness can be more accurately 
accounted for by using the offset correction factors. 

All the results presented in this paper are for bending vibration 
in a plane normal to the plane of rotation. Frequencies for vibration 
in any other plane can be easily determined from these results, however, 
by means of a simple formula proposed by Lo and Renbarger (ref. 5) - 


CONCLUDING REMARKS 


To summarize, it has been shown that the Southwell approach, based 
on nonrotating mode shapes, provides a reasonably accurate means of pre- 
dicting frequency changes due to rotational speed for most current heli- 
copter blades. It has been found, however, that a single Southwell 
constant for each mode cannot yield accurate results for all cases, 
particularly if the mass distribution of the blades is quite different. 
In order to aid the designer, charts have been presented which permit 
the rapid estimation of nonrotating frequencies, of Southwell constants, 
and of the offset correction coefficients for Southwell constants, from 
which reasonably accurate frequencies can be readily obtained for rota- 
ting blades. In example applications, the method gave good estimations 
of the bending frequencies of actual rotor blades. 


Langley Aeronautical Laboratory, 

National A.dvisory Committee for Aeronautics, 
Langley Field, Va., June l8, 195^* 



8 


MCA EM L54G02 


REFERENCES 


1. Liebers, F.: Contribution to the Theory of Propeller Vibrations. 

NACA TM 568 , 1950 . 

2. Prewitt, R. H., and Wagner, R. A.: Frequency and Vibration Problems 

of Rotors. Jour. Aero. Sci., vol. no. 10, Aug. 19^0, pp. 444-450. 

3 . Loj Hsu, and Renbarger, Jack L.: Bending Vibrations of a Rotating 

Bean. Proc . First U. S. Nat. Cong. Appl. Mech. (Chicago, 111., 

1951 ), A.S.M.E., 1952 , pp. 75 - 79 . 



MCA EM L54G02 


9 


TABLE I 


EXACT 


AND 


ESTIMATED FREQUENCIES FOR SEVERAL 
MANUFACTURED BLADES 


m 




MODE 

RAOIANS/SEC 

wr, radians/sec 

EXACT 

EST. 

EXACT 

EST. 

1st 

17.3 

17.4 

49.2 

47.7 

2nd 

48.5 

50.0 

86.8 

85.7 

3rd 

95.5 

101.0 

137.0 

137.8 



1st 

21.6 

2U 

50.6 

49.2 

2nd 

58.9 

60.5 

92.4 

92.2 

3rd 

M2.1 

122.0 

148.0 

154.0 



1 St 

21.9 

21.1 

74.0 

78.3 

2nd 

63.7 

59.5 

132.0 

134.4 

3rd 

126.0 

125.5 

200.0 

207.5 



n 

1 St 

13.4 

14.6 

37.9 

37.8 



2nd 

43.7 

41.6 

71.0 

70.3 

.J , 

— L_ 

3rd 

94.9 

94.5 

125.0 

124.0 


ROOT TIP ROOT TIP 
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Figure 5*- Effect of rotational speed on the bending frequencies 

of a uniform cantilever beam. 
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Figiire 6.- Bending frequency coefficients for hinged beams with 

linear mass and stiffness distributions. 
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Figure 10.- Zero-offset Southwell constants Kq^ for cantilever beams 
with linear mass and stiffness distributions. 











